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Abstract 

We propose a mechanism which produces periodic variations of the degree 
of predictabihty in dynamical systems. It is shown that even in the ab- 
sence of noise when the control parameter changes periodically in time, below 
and above the threshold for the onset of chaos, stochastic resonance effects 
appears. As a result one has an alternation of chaotic and regular, i.e. pre- 
dictable, evolutions in an almost periodic way, so that the Lyapunov exponent 

is positive but some time correlations do not decay. 
05.45. -hb,05.40.-Fj 



1 



Typeset using REVTgX 



The mechanism of stochastic resonance was initially introduced as a possible explanation 
of the long time climatic changes In the last years it has been used in a wide class of 

systems in physics and biology such as analogical circuits P, ring laser 0, neurology 
bistable systems Qjl^, systems with colored noise see Ref. for a recent review. 

The phenomenon can show up in bistable systems with a periodic forcing and a random 
perturbation. A typical example is the evolution generated by the stochastic differential 
equation 

dx dV{x,t) 



dt dx 



+ V2ar] (1) 



where \^ is a time periodic double well potential 

4 2 

= - y + Axcos(cjt), (2) 

and T] white noise. 

It can be shown that there exists a range of values of A, a and u where the jumps 
between the two oscillating wells are strongly synchronized, as a consequence of a sort of 
resonance between periodic forcing and random perturbation. 

Systems showing stochastic resonance are, in some sense, intermediate between regular 
and irregular ones, since they are described by a random process - the jumps do not follow a 
deterministic rule - which, nevertheless, exhibits a certain degree of regularity. For instance, 
the X time-correlation does not decay. 

This letter shows that a similar behavior can arise in deterministic systems close to the 
onset of chaos when the control parameter varies periodically in time. Under appropriate 
conditions, the time evolution shows an alternation of regular and chaotic motion strongly 
synchronized with the time variation of the control parameter. The presence of deterministic 
chaos plays the role of the random perturbation, so that it would be more correct to speak 
of 'chaotic resonance' rather then of stochastic resonance. Chaotic resonance seems to be 
present in natural phenomena, such as the time evolution of weather, which is governed by 
a set of non-linear equations which surely exhibits deterministic chaos. Nevertheless, one 



observes some elements of regularity such as high predictability during summer, at least in 
the tempered regions as the mediterranean countries, and very poor predictability during 
winter. 

Although the alternation of seasons cannot be described by low dimensional systems, 
some qualitative features can be captured by toy models, which can be useful as a first step 
toward of the comprehension of the mechanism producing periodic variations of predictabil- 
ity in short and long term climate phenomena. 



We have chosen to analyze the Lorenz model [0 which is the first geophysical dynamical 
system where deterministic chaos has been observed. We consider the original differential 
equations 

dx/dt= 10 (y- a;) 

dy/dt= -xz + R{t)x-y (3) 
dz/dt = xy — ^z 
where the control parameter has a periodic time variation: 

R(t) = Ro-A cos(27rt/T). (4) 

The Lorenz model describes the convection of a fluid heated from below between two layers 
whose temperature difference is proportional to the Rayleigh number R. In our case, the 
periodic variations of R roughly mimic the seasonal changing on the solar heat inputs. 

In order to get stochastic resonance effects without noise, the average Rayleigh number 
Ro is assumed to be close to the threshold Rcr = 24.74 for the transition from stable fixed 
points to a chaotic attractor in the standard Lorenz model. The value of the amplitude A 
of the periodic forcing should be such that R{t) oscillates below and above Rcr- For very 
large T, a good approximation of the solution is given by 



x{t) = y{t) = ±y-(i?(t) - 1) z{t) = R{t) - 1 (5) 

which is obtained by the fixed points of the standard Lorenz model by replacing R by R{t). 
The stability of this solution is a rather complicated issue, which depends on the values of 



Rq, a , and T. For instance, when Rq = 23.3 and A = 4, we numerically found that the 
solution is stable for any value of T, although R{t) can become larger than Rcr- 

On the other hand, it is natural to expect that if Rq is larger than Rcr the solution is 
unstable. In this case, for A large enough (at least Rq — A < Rcr) one observes a mechanism 
similar to that of the stochastic resonance in bistable systems with random forcing. As in 
the case of the stochastic resonance we have a periodic variation in the dynamics (the control 
parameter) and the chaos plays the role of the noise. The value of T is crucial: for large T 
the systems behaves as follows. It is convenient to call 

T„ ~ nT/2 - r/4 (6) 

the times at which R{t) = Rcr- For < t < Ti, the control parameter R{t) is smaller 
than Rcr so that the system is stable and the trajectory is close to one of the two solutions 
(I). For Ti < t < T2, one has R(t) > Rcr and both solutions (H) are unstable so that the 
trajectory in a short time relaxes toward a sort of 'adiabatic' chaotic attractor. The chaotic 
attractor smoothly changes at varying R above the threshold Rcr, but if T is large enough, 
this dependence can be neglected in a first approximation. However, when R{t) becomes 
again smaller than Rcr, the 'adiabatic' attractor disappears and, in general, the system is 
far from the stable solutions (^. But, since they are attracting, the system relaxes toward 
them. If the half-period Ti is much larger than the relaxation time tc, in general the system 
follows one of the two regular solutions (|^) for T2n+i < t < T2„+2- However, there is a small 
but non-zero probability that the system has no enough time to relax to (^) and that its 
evolution remains chaotic. Figure |l] shows the time evolution for T = 300 (a) and T = 1600 
(b). They provide a unambiguous numerical evidence that the jumps from the chaotic to 
the regular behavior (and the contrary) are well synchronized with R{t), with probability 
close to 1 when the forcing period T is very long, as in Fig. lb. On the other hand, for small 
value of T the system often does not perform the transition from the chaotic to the regular 
behavior, see Fig. la. 

It is worth stressing that the system is chaotic. In both cases, in fact we found numerically 
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that the first Lyapunov exponent is positive, although the correlation function of the variable 
z does not decay. This is due to the presence of strong correlation between the regular 
intervals. 

Figure shows the probability distribution of the lengths of the irregular interval. One 
observes peaks around T/2, 3T/2, 5T/2 • ■ while the envelope of the probability distribution 
decreases exponentially. This feature can be easily explained. At t = {n = 1, 2 ■ ■ ■) the 
system will be in some part of the 'adiabatic' chaotic attractor. The phase space is divided 
into two regions Qi and Q2 such that if x(T2„) is contained in Qi the trajectory during the 
following half-period will be very close to one of the two solutions @. On the other hand 
the points x(T2„) contained in Q2 generate trajectories which remain far from (j^). Calling 
vr the measure of the region Q2 and noting that in the irregular intervals the correlations 
decay very fast, it follows that the probability, that the lengths of the irregular interval 
is close to T2n+i is Pn — = exp(— en) with c = — Inp. 

This feature has been observed in many other systems exhibiting stochastic resonance 

The probability of jumping a regular interval, vr, decreases with the period of the forcing 
T, of course. Figure ^ shows that in the Lorenz model (^), the probability P{T) to have an 
irregular interval longer than T decreases as: 

POO 

P[T) = p{t) dr ~ e-"^ (7) 

where p(r) is the probability distribution of the length of the irregular interval. 

Without entering in the details, we briefly discuss the effect of a random forcing, of 
strength a, in the case where R(t) — Rcr changes sign during the time evolution but the 
solutions (!]), in the absence of the noise, are stable. In practice, we consider the Langevin 
equation 

dx/dt = 10 {y — x) + \^2ar]i 
< dy/dt = -XZ + R{t) x-y + V2ar]2 (8) 
dz/dt= X y — ^ z + ^/2o- 1]^ 

5 



where rii{t) are uncorrelated white noises i.e. < Tji{t)rij{t') >= SijS(t — t'). 

The numerical study of the model reveals a phenomenology very close to the original 
stochastic resonance For small values of cr one has the same qualitative behavior 

obtained at cr = 0, while for cr slightly larger than a critical value acr one has an alternation 
of regular and irregular motions. Now the Lyapunov exponent, computed treating the noise 
as an usual time-dependent term, is negative, i.e. two trajectories, initially close, with the 
same realization of the random forcing do not separate but stick exponentially fast. We 
stress that the Lyapunov exponent computed in the above method is neither unique nor the 



most physically relevant characterization of the complexity of noisy systems |jT7|. 

It is not difficult to give a rough argument for the above features. In the time interval 
where R{t) < Rcr, because of the random noise, the distance 6 between the state of the 
system x and the solutions (^) is 0{^/o^). During the half-period T2n+i < t < T2„+2, the 
typical distance 5 grows exponentially: 

5{t) ~ v^e^(*-^^"+i) (9) 

Very roughly, c is related to the largest real part of the eigenvalues of the stability matrix 
computed along the solutions (^. Calling L the size of the 'adiabatic attr actor', if the 
strength of the random forcing is large enough, i.e. 

a >(7e, ~ L'^e-^^ (10) 

the system can jump into the 'adiabatic attractor' at a time between T2„+i and T2n+2 and 
one has the same behavior shown in Fig. lb. 

This feature is quite similar to the original stochastic resonance, as the central role is 
played by the forcing term. Let us stress that the critical value acr decreases very quickly 
with the period T. 

In conclusion we have shown that the phenomenology of the stochastic resonance can 
appear in a dynamical system even in the absence of a random perturbation, when there 
is a periodic time variation of the control parameter around the onset of chaos. Instead of 



the two minima in the double well potential considered by the original stochastic resonance, 
one has two dynamical states of the system: chaotic and regular. The role of the noise is 
played by the chaotic evolution itself. It is worth noting that one needs that the period T 
of the control parameter variations should be much larger than the internal relaxation time 
tc toward the regular solution of the unperturbed system. 

Stochastic resonance in chaotic systems has relevant consequence on the predictability 
problem. It shows that the predictability time is not trivially related to the Lyapunov 
exponent if T is large enough. During the regular intervals, one has an almost perfect 
predictability while in the irregular intervals the predictability time is given by the inverse 
of the Lyapunov exponent. Moreover, we have shown that there is a non-zero probability 
(vanishing when T ^ oo) to skip a regular interval. Using a pictorial language, we could 
say that the regular interval corresponds to the summer evolution, while the irregular one to 
winter. Although the Lorenz model is too naive for any attempt of a realistic description, it 
allows us to reproduce some important features of weather forecasting which motivated our 
work: the forecasting is limited up to a time proportional to the inverse Lyapunov exponent 
of the system during winter; there is a very high predictability in summer; there is a small 
but not negligible probability to have very bad summers (jumps of the regular intervals) 
where the weather is unpredictable. 
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FIGURES 

FIG. 1. Model with A = 4., Ro = 25.5. z as a function of t/T for T = 300 (la) and T = 1600 
(lb). 

FIG. 2. Model with A = A, Rq = 25.5. Probability density, p, to have an irregular interval 
T = At/T, for T = 300. 

FIG. 3. Model with A = 4, Rq = 25.5. Probability of jumping one or more regular intervals, 
P, as a function of the forcing period, T. 
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